Conductance oscillation due to the geometrical resonance in FNS double junctions 
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We theoretically analyzed the Andreev reflection in ferromagnetic metal / nonmagnetic metal 
/ superconductor double junctions with special attention to the electron interference effect in the 
nonmagnetic metal layer. We showed that the conductance oscillates as a function of the bias 
voltage due to the geometrical resonance. We found that the exchange field and therefore the 
spin polarization of the ferromagnetic metal can be determined from the period of the conductance 
oscillation, which is proportional to the square-root of the exchange field. 



Recently much attention has been focused on the An- 
dreev reflection (AR) in ferromagnetic metal (FM) / su- 
perconductor (SC) contacts3 u 2 i 2 ia S i S i L2i2>lS since the spin 
polarization of conduction electrons is measured through 
the suppression of the conductance below the supercon- 
ducting gap. This method is called point contact An- 
dreev reflection (PCAR) spectroscopy. 

On the other hand, the quasiparticle (QP) interfer- 
ence in nonmagnetic metal(NM) / SC junctions has 
been extensively studied in the pas t 11 ' 12 ) 13 ' 14 ) 15 ' 16 ) 17 ' 18 . 
As shown in Refs . 12 i 13 ' 14 , the interference of QPs in 
the SC layer brings about the oscillation of the den- 
sity of states against the QP energy, which is called a 
Tomasch oscillation. The interference in the NM layer 
also brings about the oscillation of the density of states 
in the NM laye r 11 ! 15 , which is known as the deGennes- 
Saint- James bound state or the McMillan- Rowcll oscilla- 
tion. Nesher and Koren measured the dynamic resistance 
of YBa 2 Cu 3 6 . 6 / YBa2Cu2.55Feo.45Oj, / YBa 2 Cu 3 6 . 6 
junctions and determined the renormalized Fermi veloc- 
ity of QPs in the YBa2Cu2.55Feo.45Oy layer from the pe- 
riod of the McMillan- Rowell oscillation 16 . 

In this paper, we theoretically analyze the Andreev re- 
flection in a FM/NM/SC double junction system with 
special attention to the electron interference effect in the 
NM layer. Following the work of Blonder, Tinkham, and 
Klapwijk (BTK)i2, we solve the Bogoliubov-de Gennes 
(BdG) 20 equations and calculate the conductance. We 
show that the conductance due to the Andreev reflection 
oscillates as a function of the bias voltage because of the 
geometrical resonance predicted by deGennes and Saint- 
James. We obtain the analytical expression of the prob- 
ability of the Andreev reflection under the Andreev ap- 
proximation and find that the period of the conductance 
oscillation is proportional to the square-root of the ex- 
change field. Therefore, we can determine the exchange 
field and therefore the spin polarization of the FM layer 
from the period of the conductance oscillation. 

The system we consider is comprised of the 
FM/NM/SC double junctions shown in Fig. [Ha). The 
current flows along the x-axis, and the interfaces between 
FM/NM and NM/SC are located at x = and x = d, 
respectively. The system is described by the following 



BdG equatio: 



,20. 



Ho — h(x)a 
A*(x) 



A(x) 
-Ho — h(x)a 
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(1) 

where Ho = —(h 2 /2m)V 2 + V(x) — piF is the single parti- 
cle Hamiltonian, E is the QP energy measured from the 
Fermi energy fip, V(x) is the interfacial barrier—, and 
(j = +(— ) represents the up- (down-) spin band. The ex- 
change field function h(x) is given by h(x) — ho[l — ®(x)] 
where ho represents the exchange field in the FM layer 
and Q(x) is the Heaviside step function. We employed 
the two-band Stoner model for the FM layer for simplic- 
ity. The superconducting gap function is expressed as 
A (2) = A O(x — d), where A represents the supercon- 
ducting gap in the SC layer. We assume that the system 
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FIG. 1: (a) Schematic diagram of a FM/NM/SC double junc- 
tion. An NM with a thickness of d is sandwiched by FM and 
SC layers, (b) Schematic diagrams of energy vs. momen- 
tum of the FM/NM/SC double junction for a spin-up incident 
electron are shown. The open circles denote holes, the filled 
circles electrons, and the arrows point in the direction of the 
group velocity. The incident electron with up-spin is denoted 
by 0, along with the resulting scattering processes: Andreev 
reflection (1), normal reflection (2) at the FM/NM interface, 
transmission to the NM (3, 4) and reflection at the NM/SC 
interface (5, 6), transmission as a electron-like quasi-particle 
to the SC (7) and that as a hole- like quasi-particle (8). (c) 
Schematic diagrams of energy vs. momentum in the FM layer 
for a spin-down incident electron are shown. 
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has translational symmetry in the transverse (y and z) 
direction, and therefore the wave vector parallel to the 
interface k\\ = (k y ,k z ) is a conserved quantity. 

The general solutions of the BdG equation ((T|) in the 
FM (NM) layer are given by 

* ±fe + (r) = (T) e ±ik ™(™h« x S ( r|| ) ( 2 ) 



±fc 



FM(NM), 



(r) 



f! e ±ifc ™ ( NM),„- Sfc (r,,), (3) 



In the NM layer (0 < x < d), we have 



«Sf.(r) 



Ct.a- Lie " 



and in the SC layer (x > d) 



~ W 1 x I 



S fc|l (ry), 



(7) 



where Sfc., (ry) represents the eigen function in the trans- 
verse direction in the fey channel and ipL^j a is the x 
component of the wave number of an electron (hole) with 
tr-spin defined as k^ M a = y//J-F ± E + aho — E\\ 



and A: NM — — 
the SC layer, we have 



y/fiF ± E — E\\, where E\\ 



2m 



111 



«0j e±<fc + a x Sfcu(p||)) (4) 



*±^={2) e±ik " cXSk ^ (5) 

where uo and vq are the coherence factors expressed as 

, and kgQ is the x com- 
ponent of the wave number of an electron- (hole-) like QP 

2m 
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defined as kg C — 



± VE 2 - A 2 - En 



The wave function of the FM/NM/SC double junction 
is given by the linear combination of the above general 
solutions. Let us consider the scattering of an electron 
in the fey channel with c-spin injected into the NM from 
the FM, the eight processes shown in Fig[JJ(b) are active. 
Therefore, the wave function in the FM layer (x < 0) 
takes the form 



",T k I ^ I 



Sfc„ (ry). 



(6) 



llTkii Q) e -ik- sa , 



(8) 



Sfe^ry). 



The coefficients a a , kl] ,b a . kll ,c a ,k\\ > ,<x<r,k\\ $a,k\\ > Co-,fc n , 
and X<7,k|i are determined by matching the wave func- 
tion at the boundary of the contact x — and d. Fol- 
lowing the BTK theoryi^ the probabilities of the AR 
and the normal reflection are given by A^kn (E) — 

( k FM,a/ k FM,aX,k^,h and B a>h ^{E) = re- 
spectively. Since we assume that the temperature is 
zero, the conductance at bias voltage V is given by 

G = f Ea.fey i 1 + 4r,fc|| ( e ^) ~ ( eV )] > where we as- 
sume that the voltage drop occurs at the NM/SC inter- 
face for simplicity. Below the superconducting gap, i.e., 
eV < A , the probabilities A^^^E) and B aj k^{E) sat- 
isfy the relation that 1 — B a ^ (E) = A ai k» (E) and then 
we have 



(9) 



Let us first consider the most idealistic case where the 
interfacial scattering potential, V{x), is assumed to be 
zero. This simplification enables us to obtain the analyt- 
ical expression of A a ^ {E) under the Andreev approxi- 
mation. 



A*m(E) 
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(10) 



where we introduced the normalized parameters rj = 
ho/up, ( = and e = E/A . We note that Eq. 

()10|1 contains trigonometric functions in the denomina- 
tor. For the FM/SC junction, i.e., d—0, the trigono- 
metric functions become constant and Eq. (|10[) repro- 
duces the deJong and Beenakker results of the zero-bias 
conductance!. For the FM/NM/SC junctions with d ± 0, 



the trigonometric functions in Eq. (|10[) give rise to the os- 
cillation of the conductance against the bias voltage. The 
origin of the conductance oscillation is the interference of 
electrons in the NM layerii^. In the FM/NM/SC dou- 
ble junctions, the injected electron propagates across the 
NM layer to the interface as an electron (3 in FigfT]) and 
is scattered into a hole (6 in FigfT]) by the superconduct- 
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FIG. 2: (a) Probability of the AR A afi (eV) for the 
FM/NM/SC double junction with d=l/xm is plotted against 
the bias voltage V . The exact numerical results and the value 
of Eq. (fT0)l are plotted by lines and circles, respectively. The 
value of the exchange field is taken to be h = 0.3, 0.6, 0.9 /j,f 
from top to bottom, (b) Same plot for d=10^im. 



ing gap. The superconducting gap can pair an excited 
electron with an electron inside the Fermi sea, leaving 
a hole excitation. The hole propagates back across the 
NM layer; however, it cannot interfere with the original 
electron. In order for interference to occur the hole must 
be reflected at the FM/NM interface, propagate to the 
NM/SC interface, be scattered into the electron state (5 
in FigfJ} by the superconducting gap, and propagate as 
an electron in the NM layer. It can interfere with the 
original electron (3 in Fig[T]). This interference produces 
an oscillation of the conductance against the bias volt- 
age, and the period of the oscillation is determined by 
the thickness of the NM layer. 

In order to analyze the interference effect on the con- 
ductance oscillation, we consider the AR probability 
-^-o-.feii {E) of the one-dimensional system; i.e., only the 
transverse channel with fen = is considered. In Fig. 
[2] (a), H (b) we plot the probability (eV) of the 

FM/NM/SC double junctions with d=l/im and 10/zm, 
respectively, as a function of the bias voltage. Since 
Eq. (|10p is an even function of normalized value of the 
exchange field rj, A^^^ {eV) is independent of the spin di- 
rection tr; i.e., A^^ (eV) = Aufe.. (eV). The Fermi energy 
and the superconducting gap are assumed to be \ip = 3.8 
eV (kp — l.oA ) and A = 1.5meV, respectively. The 
exact numerical results and the approximate values of 
Eq. (fTU|) are plotted by lines and circles, respectively. 
The value of the exchange field is taken to be ho = 0.3, 
0.6, 0.9 /xp from top to bottom. 

As shown in Figs. M (a) and H (b), Eq. flU]) and 
therefore the Andreev approximation are valid for all 
values of the exchange field. According to Eq. (flU)) . 
the period of the oscillation is determined by the condi- 
tion that k + — k~ = rnr, where n is an integer. Since 
k ~ fcp(l ± E/2/ip), the period is obtained as 



AVi 
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2m e d 



(11) 



which is inversely proportional to the thickness of the 
NM layer, d. For the one-dimensional system, the period 
AVid is independent of the exchange field of the FM layer 



as shown in Figs. [2] (a) and [2] (b), and in Eq. lfTTj) . How- 
ever, as we shall show later, the period of the conductance 
oscillation due to the geometrical resonance depends on 
the exchange field of the FM layer because the number 
of fe|| channels available for the AR is restricted by the 
exchange field. 

The period of the oscillation of A^k* i e V) with finite fey 
is given by AVfe,, ~ J^ ed ~ E \\ ■ From E( h© the 
conductance is obtained by summing up A^k.. (eV) for all 
available fey . Since the spin of the Andreev reflected hole 
is opposite to that of the incident electron, the maximum 
value of fcy and therefore E\\ is limited by the exchange 
field ho as max-Ey = fip — h . We assume that oscilla- 
tions of A a .k^ (eV) with different periods cancel out each 
other and the period of the sum fe A a _k }l ( e V) is de- 
termined by the shortest period. Thus, the period of the 
conductance oscillation of the three-dimensional system 
is obtained as 



AV 3B ~ minAFfc.. 



hn 



(12) 



In Figs. [3](a),[3](b), and[3](c) we plot the conductance 
of the FM/NM/SC junction, Gfns, normalized by that of 
the FM/NM/NM junction, G F nn, against the bias volt- 
age. As shown in Fig. 0(a), the oscillation due to the ge- 
ometrical resonance does not appear in the conductance- 
voltage curve if the thickness of the NM layer, d, is less 
than or of the order of nm. The conductance-voltage 
curve is indistinguishable from that of the FM/SC junc- 
tion. Hence, we can use the conventional PCAR analysis 
for a FM film, the surface of which is coated by a thin 



(a) d = 1nm 








— h 


=0 




h 


=0.3hf 




h 


=0.6 M f 




h 


=0.9 M f 


j 









C 


1 2 
V[mV] 


3 




(c) d = 10|im 




2 




h„=0 

h =0.3HF 

— h -0.6^F 


z 
z 




- - h -0.9^F 


a 

- 1 






z 






a 












"0 1 2 
V [mV] 


3 



(b) d = 1|im 






h =0 

h -0.3|j F 

h =0.6MF 

h =0.9^ F 











E 



0.01 L 

0.01 



(d) d = 10(im 



0.1 

ho M 



FIG. 3: (a) The conductance of the FM/NM/SC double junc- 
tion, Gfns with d=lnm is plotted against the bias voltage. 
The conductance is normalized by that of the FM/NM/NM 
junction, Gfns- (b) Same plot for d—lfim. (c) Same plot for 
d=10fim. (d) The period of the conductance oscillation of the 
FM/NM/SC double junction with d=10^m is plotted against 
the exchange field ho. 
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FIG. 4: (a) The normalized conductance, Gfns/Gfnn for 
d=lfim, Z—0.2 is plotted against the bias voltage, (b) Same 
plot for d=10^m, Z=0.2. 

(less than a few nm) NM layer. 

Figures [3] (b) and [3] (c) show the conductance- voltage 
curves for the FM/NM/SC double junctions with d = 
1/xm and d = 10/im, respectively. One can see that the 
period does depend on the exchange field, ho, in the FM 
layer as well as the thickness of the NM layer. The pe- 
riod is a decreasing function of the exchange field. We 
can easily confirm that the period is proportional to the 
square-root dependence of the exchange field by looking 
at Fig. [3] (d). The exact numerical results (filled cir- 
cles) agree well with Eq. (fT2"]) (dotted line). The results 
suggest that we can determine the exchange field and 
therefore the spin polarization of the FM layer from the 
period of the conductance oscillation. 

Next we consider the effect of the interfacial scattering 
potential at the NM/SC interface. We assume that the 
interfacial potential is represented by the delta-function 
as V(x) = {h 2 kFZ/m)8{x~d), where Z is the dimension 
less parameter which characterize the strength of the in- 
terfacial scattering potential. For simplicity, we neglect 
resistances of the FM/NM interface and the FM layer 
which do not change the oscillation period of the conduc- 



tance but reduce the amplitude of it. In Figs[4](a) and[4] 
(b), we show the normalized conductance- voltage curves 
for junctions with d=l/zm and d=10/im. The parameter 
Z is assumed to be 0.2&. Because the conductance os- 
cillation is due to the geometrical resonance in the NM 
layer, the period of the oscillation is not affected by the 
interfacial scattering potential and is given by Eq. lfT!?]) . 

In the present analysis we employed the simplest BdG 
approach and consider the clean FM/NM/SC junctions 
with perfect interfaces. In the real experiments the con- 
ductance oscillation we predicted might be smeared out 
due to the interface roughness and imperfections. How- 
ever, recent advances in fabrication technology enables 
us to fabricate epitaxial FM/NM/SC trilayers of high 
quality^. We expect that the conductance oscillation 
we predicted can be observed in such epitaxial trilay- 
ers. For further understanding of the transport prop- 
erties of FM/NM/SC trilayers, we have to take into ac- 
count the effects of finite mean free path, band structures, 
and selfconsistent determination of the electron's distri- 
bution function and electric-pottential, which is beyond 
the scope of this Brief Report. 

In summary, we studied the conductance oscillation 
due to the geometrical resonance in a FM/NM/SC double 
junction theoretically. We showed that the conductance 
due to the Andreev reflection oscillates as a function of 
the bias voltage due to the geometrical resonance. We 
found that the exchange field and therefore the spin po- 
larization of the FM layer can be determined from the pe- 
riod of the conductance oscillation because the period of 
the conductance oscillation is proportional to the square- 
root of the exchange field. 

The authors would like to acknowledge the valuable 
discussions they had with K. Matsushita and N. Yokoshi. 
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